Can disorder alone destroy the e' g hole pockets of Na^CoC^? 
A Wannier function based first-principles method for disordered systems 
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We investigate from first principles the proposed destruction of the controversial e' g pockets in the 
Fermi surface of Na I Co02 due to Na disorder, by calculating its fc-dependent configuration-averaged 
spectral function {A(k,oj)). To this end, a Wannier function-based method is developed that treats 
the effects of disorder beyond the mean field. Remarkable spectral broadenings of order ~leV are 
found for the oxygen orbitals, possibly explaining their absence in the experiments. In contradiction 
with the current claim, however, the e' g pockets remain almost perfectly coherent. The developed 
method is expected to also generate exciting opportunities in the study of the countless functional 
materials that owe their important electronic properties to disordered dopants. 



PACS numbers: 74.70.-b, 71.15.-m, 71.18.+y, 71.23.-k 

Like most transition metal oxides, sodium cobaltates 
obtain their important properties via the introduction of 
dopants (Na) between the layered structures of oxygen 
and transition metal atoms. Around x = 0.3, Na^OoC^ 
develops unconventional superconductivity under hydra- 
tion pQ, with evidence of a nodal order parameter [2]. 
From x < 0.5 to x > 0.5 it changes from a paramag- 
netic metal to a Curie- Weiss metal, while it is a charge 
ordered insulator at x = 0.5 [3J- At high doping, the 
combination of high thermopower and high conductiv- 
ity is observed [4 , together with A-type antiferromag- 
netism 5J. In addition, various Na orderings have been 
observed throughout the entire phase diagram [BJ. This 
rich variety of behaviors has thus attracted intensive re- 
search activity. 

Nevertheless, even the most basic starting point for an 
understanding is still under serious debate, namely the 
low-energy electronic structure near the chemical poten- 
tial that controls most of the above mentioned remark- 
able properties. Early density functional (DFT) calcula- 
tions 7 within the local density approximation (LDA) 
predicted the existence of a central a g hole pocket, sur- 
rounded by 6 e' g hole pockets. Angle resolved photoemis- 
sion spectroscopy (ARPES) experiments [8, 9J measured 
the central a g pocket, but found the e' g bands to be below 
the Fermi surface. Shubnikov-de Haas measurements [TU] 
observed two pockets, but the assignment to e' g was con- 
cluded incompatible with the specific heat data pj] . On 
the other hand, the presence of the second type of pocket 
was reconfirmed by Compton scattering |12j and assigned 
to e' g . 

Despite the controversial status on the experimental 
front, various theoretical efforts have been made to in- 
vestigate the alleged absence of the e' g pocket. Surface 
effects were suggested to suppress the e' g pocket [13] un- 
der hydroxyl contamination. The other usual suspect of 
many-body correlation was investigated within dynami- 



cal mean-field theory (DMFT) by several groups. How- 
ever, the results were inconclusive as the e' g pockets were 
found to either grow or shrink depending sensitively on 
the crystal field splitting [T3j. Some researchers [TS] ar- 
gued that the e' g pockets should not exist according to 
the specific heat data, while others [16] concluded from 
an extensive study of the crystal field that the e' g pockets 
cannot be removed via local correlation. 

Recently, an intriguing alternative resolution was pro- 
posed 17J. It was argued that the random positioning 
of the Na intercalants alone can introduce strong disor- 
der effects that mask the e' g pockets from the ARPES 
experiments. This physically plausible picture, if proven, 
would not only enable a new resolution to reconcile the 
various theoretical and experimental observations, but 
would also introduce important novel physics missing in 
current considerations. Unfortunately, verification of this 
intriguing proposal presents a great challenge to the cur- 
rent capability of the first-principles theories. To un- 
ambiguously resolve the complete disorder effects with a 
good fc-space resolution, the disorder-induced self-energy 
cannot be assumed local a priori. Especially with the 
potential localizationpj] near the pockets, the spatial 
fluctuations need to be taken into account, beyond the 
current state-of-the-art mean-field treatments [T5H2D] . 

In this Letter, we examine the proposal of disorder- 
induced destruction of the e' g pockets in Na^OoC^, us- 
ing x = 0.3 as a representative case, by developing a 
first-principles Wannier function-based method for the 
evaluation of electronic structure of disordered materials. 
Remarkable spectral broadenings (~leV) of the oxygen 
bands arc found that provide a natural explanation for 
the missing oxygen bands in ARPES studies. However, 
in contradiction with the current claim [17] , we only find 
a negligible influence of disorder on the e' g pockets, in- 
capable of masking them from the ARPES experiments. 
The inclusion of nonlocal disorder-induced self-energy in 
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our method is expected to generate exciting opportuni- 
ties in the studies of countless modern functional ma- 
terials, including doped transition metal oxides, dilute 
magnetic semiconductors, and intercalated graphites, to 
name but a few. 

The electronic band structure of a disordered system is 
determined by the configuration-averaged spectral func- 
tion: (A(k, w))=X^ c PcA c (k, ui) of momentum k and en- 
ergy u, in which configuration c is weighted by its prob- 
ability P c . By treating the disordered configurations 
within the supercell approximation, the spectral func- 
tions A c (k, of) can be calculated directly from the su- 
percell eigenvalues and eigenvectors, by applying the un- 
folding method of [21] . While conceptually clean, the ap- 
proach is in practice too computationally expensive for 
most applications due to the large size of the supercells 
required to incorporate the characteristic length scale of 
localization and impurity scattering. The goal of this 
work, therefore, is to drastically reduce the computa- 
tional cost by developing an effective low-energy disor- 
dered Hamiltonian using first-principles Wannier func- 
tions. 

The Hamiltonian of an arbitrary configuration of N im- 
purities, positioned at (x%, xn), can be exactly rewrit- 
ten as 
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where H° denotes the Hamiltonian of the system with no 
impurities, A^^—H^^-H , denotes the linear influence 
of the impurity at xt and A^'^'^iJ^'^-A^-A^')- 
H° denotes the two-body correction of a pair of impuri- 
ties at (xi,Xj), etc. Naturally, low-order Greens functions 
like the spectral functions here are not too sensitive to 
higher order corrections and the computation can be sim- 
plified accordingly. In practice we found that it is already 
highly accurate to keep only the linear influence of the 
impurities, as demonstrated below. It is important to 
note that a cutoff here by no means implies a truncation 
in the electronic multiple scattering processes among the 
impurities. Also, while Eq. ([I]) resembles the well es- 
tablished cluster expansion method [55] of the formation 
energy of large-sized impurity clusters, it encapsulates 
the influence of disorder in the entire low-energy effec- 
tive Hamiltonian. 

The construction of the effective Hamiltonian explic- 
itly consists of three steps. First, DFT calculations are 
performed for the undoped normal cell, and for each type 
of impurity in a large supercell containing one impurity. 
Second, the influence of each impurity (located at x m in 
the first normal cell) is extracted, A^" 1 ) -H° , in the 

basis of Wannier functions. Similar to the induced forces 
in the typical frozen phonon calculations, the partition 
of the influence of the impurity from its superimages is 
necessary [53]. Third, the effective low-energy Hamilto- 
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FIG. 1: (color online) Isosurface plot |(a;|m)| = 
0.09(bohr) -3 ' 2 of selected Wannier functions, colored accord- 
ing to the phase of the complex functions. 



nian corresponding to each disorder configuration with 
N impurities is then assembled: 
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where r and n denote the lattice vector and the orbital 
index of the Wannier functions, and r.j and rrii denote 
the lattice vector and the type of impurity i, located at 
Xi = Ti + x mi . The notation here also implies possible 
addition of impurity orbitals and removal of orbitals due 
to vacancies. Obviously, a similar procedure can be per- 
formed if higher order contributions are desired. 

Note that only when the supercell and normal cell 
Hamiltonians are represented on the same Wannier basis 
can one legally add and subtract them. Therefore, the 
Wannier functions should be constructed from a large en- 
ergy window such that the Hilbert space is as complete as 
possible. Furthermore, the projected Wannier function 
method [24, 25 naturally ensures maximum consistency 
of the normal cell and supercell [53], and is thus more 
suitable than the maximally localized Wannier function 
method [26] which risks defining the gauge differently in 
the supercell in favor of better localization near the im- 
purity. 

For the case of Naa;Co02, we consider both possible 
locations of Na between the oxygen layers discussed pre- 
viously in the literature [17]: above Co [Na(l)] or above 
a hole in the Co sheet [Na(2)]. The low-energy Hilbert 
space is taken within [-7,3] eV consisting of symmetry- 
respecting complex Wannier orbitals of Co-<i and O-p 
characters (22/normal cell) [53], as illustrated in Fig. [I] 
The effects of impurities are extracted from three DFT 
calculations [23]: the undoped C02O4 in the normal cell 
and Na(l)Co 8 Oi 6 and Na(2)Co 8 Oi 6 in 2 x 2 x 1 supercells 
corresponding to 1 = |. Since effects of lattice relax- 
ation on the e' g pockets are found to be negligibly small 
in comparison [53], we proceed with the experimentally 
averaged structure of Nao.3Co02 in the following analy- 
sis. 

The quality of our effective Hamiltonian can be ver- 
ified by benchmarking the spectral function of a test 
system, against the standard DFT. As a highly nontriv- 
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FIG. 2: (color online) Spectral function of 

Na(2)Na(l) 2 Co 6 Oi2 obtained from (a) the full DFT 
calculation, on a basis of 2019 LAPW's (b) a single diag- 
onalization of the effective Hamiltonian, on a basis of 66 
Wannier functions, (c) Na positions in the supercell, in 
relation to triangular Co sheet. Light(dark) circles denote 
Na above(below) the Co sheet. 



ial test case we take the periodic Na(2)Na(l)2Cc-60i2(cf. 
Figj^c)), which requires a strongly "incommensurate" 
extrapolation in the partitioning (from 2 x 2 x 1 to 
\/3 x v3 x 1 cell) and the linearity (from x = ^ to 
x = 4 + |) of the influence of impurities. As shown 
in Fig. [2j our effective Hamiltonian manages to repro- 
duce the spectral function of the full DFT calculation 
with high accuracy (in particular the details around the 
gap opening), but with only a negligible fraction of the 
computational effort. (The full DFT calculation involved 
^20 self-consistent cycles on a basis of 2019 linear aug- 
mented plane waves, while our effective Hamiltonian re- 
quires only a single diagonalization on a basis of 66 Wan- 
nier functions.) Additional benchmarks exploring the po- 
tential limitations of our approximations for large extrap- 
olations are given in Ref. |23j . 

Having an accurate and efficient method to assemble 
the effective Hamiltonian of any configuration, we pro- 
ceed to evaluate the configuration-averaged spectral func- 
tions for the case of disordered Nao.3Co02. Following 
the considerations laid out in Ref. [T7], all the random 
configurations of Na are assumed comparable in proba- 
bility, except the high energy case containing two nearest 
neighboring Na atoms located at Na(l) and Na(2) sites 
[Fig. [3}d] , which is disregarded due to its low probabil- 
ity. Figure [3] shows the resulting spectral function con- 
verged with respect to the number of configurations (50) 
and their average size (~80 normal cells corresponding 
to 80 x 22 = 1760 Wannier functions). Note that it is 
necessary to consider supercells [e.g., Fig. |3][a)] of differ- 
ent sizes, orientations, and shapes in order to remove the 
effects of artificial zone boundaries of the supercell. 

A remarkable broadening of oxygen bands can be ob- 
served in Figs. [3^d) (g) , indicating a short lifetime and 
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FIG. 3: (color online) Configuration-averaged spectral func- 
tion of Nao.3Co02, showing the e' g states to be the least in- 
fluenced by disorder, (a) One of the 50 large-sized supercells 
used for configurational average, (b) A high energy configura- 
tion with Na(l) too close to Na(2). (A(k,uj)) (c) around the 
Fermi surface (d) in the full low-energy Hilbert space, where 
the bars correspond to the energy distribution curves (e) at 
fc =2rK/3, (f) at fc =r , (g) at k =F. The Wannier orbital 
character is labeled according to the color scheme of Fig. [2] 



mean free path of the quasiparticles due to strong scat- 
tering against the disordered Na atoms. This is under- 
standable considering that the Na atoms are located in 
the oxygen cages defined by the two oxygen layers, and 
thus have the largest impact on the oxygen orbitals. In- 
terestingly, this huge spectral broadening and low inten- 
sity might explain why some of the oxygen bands are not 
observed in the ARPES measurements [HI , not only 
in Na x Co02, but also in most doped layered transition 
metal oxides, where the dopants are introduced between 
the oxygen layers. 

In great contrast, much weaker effects of disorder are 
found on the Co-d orbitals [cf. Figs. [3jc)(e)(f)]. Specif- 
ically, the a g orbital picks up some k- dependent broad- 
ening near the T point, while the e' g orbitals are almost 
unaffected by the disorder. This somewhat surprising 
result probably reflects the strong screening of oxygen 
that shields the cobalt valence orbitals from the influ- 
ence of the disordered Na, and the metallic nature of 
the doped electrons that spread throughout the whole 
system. Clearly, the localization and broadening of the 
e' g and a g bands are not strong enough to shift the e' g 
pockets below the Fermi level. Unlike the a g orbital that 
points directly toward the doped Na atom (cf.: Fig. [I]), 
the e' g orbitals neither point toward the Na atom nor the 
most affected 0-p z orbitals near the Na, making them 
the least sensitive to the presence of Na intercalants. 
This could also explain the negligible effects of lattice 
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FIG. 4: (color online) Configuration-averaged spectral func- 
tion of Nao.7Co02, showing short-range order suppressing the 
spectral broadening of a g . (a) Small islands of homogeneous 
Na(l) and Na(2). (b) (A(k,uj)) around the Fermi-surface, 
(c) energy distribution curve at fco=F. The Wannier orbital 
character is labeled according to the color scheme of Fig. [2] 



relaxation around Na atoms on the e' pockets 
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viously, the disorder alone does not destroy the six e' g 
Fermi pockets of Nao.3Co02, in contradiction with the 
current claim [17] . 

On the other hand, the Na impurities do introduce an 
important physical effect on the Co-d shell, namely on 
the crystal field splitting (the relative on-site-energy) of 
the e' g and a g orbitals. Indeed, evaluated from H° and 
A*"' 1 ', the crystal field splitting is found to change from 
27 meV to -38 meV for Co atoms right below or above the 
Na intercalants, consistent with the trend estimated pre- 
viously from the density of states [2U]. Considering the 
tendency of strong orbital polarization of the many-body 
exchange interactions as demonstrated from the previ- 
ous DMFT studies [HI [20] , the combination of disorder 
and strong exchange interactions is very likely to give 
stronger scattering for the e' g orbitals. Exactly whether 
this would lead to a resolution of the highly controversial 
status of the e' g pockets of Na2,Co02, remains an inter- 
esting challenge to the theorists. 

Interestingly, our results also demonstrate a signifi- 
cant nonlocality of the disorder-induced self-energy. In- 
deed, a strongly k-dependent spectral broadening can be 
clearly observed in Fig. [3] that correlates well with the 
inverse of the band velocity rather than the energy. For 
example, at 25 meV the almost purely a g bands have 
large (~200meV) spectral broadening near the T point, 
but negligible width at k ~ |AT. Such a strong k- 
dependence of the spectral width reflects the intrinsic 
nonlocality of the self-energy, and highlights the advance- 
ment of our method over standard mean-field theories in 
which the self-energy is assumed local. 

Amazingly, related to the nonlocal self-energy, we also 
found significant effects of short-range ordering of Na im- 
purities known to be important for the Naa;Co02 [6[l22j. 
As demonstrated in Fig. |4^a) for Nao.7Co02, the exclu- 
sion of nearest neighbor Na positioning introduces auto- 
matically a strong short-range ordering of the Na impu- 
rities, due to lack of available locations at high doping 
(see also [23]). In turn, the resulting spectral function of 
the a g orbitals, for example, demonstrates stronger co- 
herence and longer lifetime, as can be seen by comparing 
Figs.gb) and Etc) with Figs. H^c) and gf). (The e' g 



orbitals near the pockets remain perfectly coherent also 
at this doping level.) This result, while physically intu- 
itive, is actually quite nontrivial, since in the mean-field 
theories the degree of disorder at 0.7 doping level should 
be exactly the same as that at 0.3 doping level. 

In conclusion, we have investigated from first principles 
the proposed destruction of the controversial e' g pockets 
in the Fermi surface of Naa;Co02 due to Na disorder. To 
this end, a Wannier function-based method is developed 
that incorporates the spatial distributions of impurities 
beyond the mean field. The new method is benchmarked 
against the full DFT calculation and shown to be efficient 
and highly accurate. Remarkable /c-dependent broaden- 
ings of the spectral function are found in the oxygen or- 
bitals due to their vicinity to the Na intercalants. How- 
ever, the effects of disorder are found to be negligible on 
the e' g orbitals. Thus, the disorder alone does not de- 
stroy the e' g pockets, in contradiction with the current 
claim [17] . Interestingly, against the mean- field perspec- 
tive, enhanced coherence is found at higher doping where 
short-range order grows stronger. Our new method is ex- 
pected to find a wide range of applications in the studies 
of the countless functional materials that owe their im- 
portant electronic properties to disordered dopants. 
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